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were first described in the context of massive 1+1-dimensional integrable quantum field theories 
and their correlation functions are related to the bi-partite entanglement entropy. We find 
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evaluating the conformal dimension of the corresponding twist field in the underlying conformal 
field theory. We find that solutions to the form factor equations are not unique so that various 
techniques need to be used to identify those corresponding to the branch point twist field we 
are interested in. The models for which we carry out our study are characterized by staircase 
patterns of various physical quantities as functions of the energy scale. As the latter is varied, 
the /3-function associated to these theories comes close to vanishing at several points between 
the deep infrared and deep ultraviolet regimes. In other words, renormalisation group flows 
approach the vicinity of various critical points before ultimately reaching the ultraviolet fixed 
point. This feature provides an optimal way of checking the consistency of higher particle form 
factor solutions, as the changes on the conformal dimension of the twist field at various energy 
scales can only be accounted for by considering higher particle form factor contributions to the 
expansion of certain correlation functions. 
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1 Introduction 



Entanglement is the most distinct and bizarre of all quantum phenomena. The idea that the 
quantum states of objects so wide apart from each other that they are not even causally con- 
nected can be entangled seems counterintuitive and has been hotly debated pQ. Although the 
interpretation of quantum mechanics and quantum entanglement continues to be discussed, the 
reality of entanglement as a physical phenomenon was finally established through the exper- 
iments of Alain Aspect and his collaborators [2|. In addition, in recent years, entanglement 
has started to reveal its powerful practical applications, specially in the context of quantum 
computation, quantum cryptology and quantum teleportation (see e.g. [3j 2]). 

Prom a theoretical point of view entanglement is intimately linked to the structure of quan- 
tum states so that developing methods to "measure" entanglement is a very efficient way to 
learn more about the fundamental properties of quantum systems. In that sense, many different 
theoretical measures of entanglement have been proposed in the literature [5j El [8j [9] of which 
the entanglement entropy [5] is just an example. In recent years much work has been carried 
out to compute the entanglement entropy of extended quantum systems with many degrees of 
freedom, such as quantum spin chains [MElIIllIIlIHIISlIIlllITllIlIISlllO] and quantum 
field theories [2H E2 [231 EH ES] . 

Consider a quantum system, with Hilbert space % = Ha <8> Hb, in a pure state \ip). The 
bi-partite entanglement entropy Sa is the von Neumann entropy associated to the reduced 
density matrix of the subsystem A, p A , defined as 

P A = Tt Hb (\^)(^\), S A = -Tx nA (p A \og(p A )) . (1.1) 

One way of interpreting 5,4 is to understand it as a measurement of the entanglement of the 
quantum state of subsystem A when the latter is considered in isolation (ignoring the existence 
of subsystem B). 

In a series of recent works involving one of the present authors |24[ [271 [251 E5] a new approach 
to the computation of (jl.ip for 1+1-dimensional integrable quantum field theories (QFT) has 
been proposed and developed. This approach takes as starting point the "replica trick". This 
consists of replacing the theory under scrutiny by a new model consisting of n non-interacting 
copies or "replicas" of the original theory. Although this might seem to lead to an unnecessary 
complication of the problem, it does in fact simplify it by producing a new theory which possesses 
an extra symmetry under cyclic permutations of the n-copies of the model. Associated to this 
symmetry there exists a special class of twist fields T and T = which have been named 
branch point twist fields. The key result is 

Tr HA p n A ~ <T(r)t(0)>, (1.2) 

that is, the trace of the density matrix of the replica theory is proportional to the two-point 
function of twist fields, r is the size of subsystem A in the quantum field theory. Therefore, 
once the trace has been computed, the von Neumann entropy follows from the identity, 

S A = - lim ^-Tr HA p n A . (1.3) 
an 

This identity does however involve a highly non-trivial analytic continuation of the function 
Try_ A p A to positive non-integer values of n which makes the limit above extremely hard to 
compute [24 [ I27 [ [28], It is sometimes more convenient to compute another type of entanglement 
entropy known as Renyi entropy [29] which is given by 
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and whose n — > 1 limit gives once more. 

In the current work we wish to take the identity (|1.2p as our main motivation to investigate 
the form factors of the twist field T. As is well known, the correlation function in (jl.2p can be 
expanded in terms of form factors of the fields involved. Therefore knowing the form factors 
allows us in principle to extract both the von Neumann and Renyi entropies of any integrable 
QFT under consideration. 

This paper is organised as follows: In section 2 we summarise the form factor programme 
for branch point twist fields [24]. In section 3 we introduce the A-sum rule [30] and the two 
models which we want to study in detail: the roaming trajectories (RT) model and the SU (3)2- 
homogenous sine-Gordon (HSG) model. In section 4 we describe in detail the construction of 
higher particle form factors for the RT-model. We place special emphasis on the existence of 
multiple solutions to the form factor equations. We successfully apply the cluster decomposition 
property to identify those that must correspond to the branch point twist field which is related 
to the bi-partite entanglement entropy of the model. In section 5 we carry out a similar analysis 
for the SU (3)2-HSG model concentrating less on the non-uniqueness of solutions and more on 
the extra challenges posed by the more complex particle spectrum of the model. In section 6 
we provide numerical results for the ultraviolet conformal dimension of the twist field in the 
RT- and HSG- model which are fully consistent with the theoretical predictions. We present our 
conclusions in section 7. 



2 Form factors of branch point twist fields 

In this paper we will be concerned with the computation of matrix elements of the branch point 
twist field T for particular models. Once the form factors are known they may be used in the 
expansion of any correlation functions involving the twist field. For example, the two-point 
function (T(O)O(r)) can be written as 

n ( k 00 \ 

<r(r)o(o)> = y: (IT /7^K r|M (0i,---A;n) 




k 

r ^2 m M • cosh t 



(2.5) 



where O represents a generic local field of the theory. Here fT^ 1 . . . ,#fc,n) is the k- 

particle form factor of the field T (similarly for the other field) which is defined as 

if 1^-^(0!, . . . , k ; n) := <O|T(O)|0i, . . . , 0, % „ ; , (2.6) 

where |0) represents the vacuum state and . . . , #fc)Jf „ i/x are the physical "in" asymptotic 
states of massive QFT. They carry indices which are double indices of the form 

Hi := (oi,Ci), (2.7) 

where on labels the particle species and Cj labels the copy number within the replica theory. 
The mass of the corresponding particle is denoted by m^. and its energy and momentum are 
parameterized by the real parameter 0i, called the rapidity. The form factors also depend on 
n, the number of replicas of the model whose entropy we want to investigate. In particular, 
for n = 1 the twist field should reduce to the identity field so that all form factors, except the 
vacuum expectation value should vanish. 
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For most QFTs the form factors defined above are only accessible perturbatively. However, 
one of the remarkable advantages of integrable QFTs is that integrability is constraining enough 
as to almost fix the functions (|2.6p completely. This feature was realized several decades ago 
|3H \32\ [33] when it was shown that form factors such as (|2.6p may be systematically obtained as 
the solutions to a set of consistency equations which only require the knowledge of the scattering 
(S) matrix of the theory under consideration as input. This solution process is sometimes 
referred to as the form factor programme. This fact has triggered an enormous amount of work 
in computing form factors in a multitude of models of integrable QFT (see |34[ 131)1 ES] for 
reviews) . 

The special nature of the twist fields T and T has however made it necessary to rethink 
the form factor programme described above from first principles in order to adapt it to replica 
theories [24] . In order to make expressions more compact, we will be dropping the explicit 
n-dependence of the form factors defined by (|2.6p . If we consider integrable QFTs without 
backscattering or bound states (theories with backscattering and bound states have been con- 
sidered in [28] and [21], respectively), then the new form factor equations are given by 

if' - M ^ +1 -(- • • , Oi, Oi+i, ■■■) = S^l +i (6 ll+1 )F^ + ^-(. . . , 9 t+1 ,9 u . . .), (2.8) 
Fl^-^{6 l + 2m,...,e k ) = F fe r|M ^2,... A,0i), (2-9) 

where 9ij = 6{ — 9j and fi i = («j, q + 1). The function Sj£ji i+1 (9a + i) is the two particle S-matrix 
of the replica theory, defined as 

S&1 (9) = ( !? Cl * Cl+1 (2.10) 

ma+iV J \ S aiai+1 (9), iff Ci = (H +1 y ' 

and S aiai+1 (9) represents the S-matrix of the original (non-replica) model. Equations (|2.8p is 
the same for form factors of other local fields, whereas equation (|2.9p is slightly different, in that 
the index fi i is involved. This difference is responsible for the fact that, on top of the standard 
kinematic poles at 9ij = in, the form factors have extra poles at 9ij = m(2n — 1) in the extended 
physical strip Q(9ij) € [0,27m). Thus we have two kinematic residue equations, which provide 
expressions for the residues at those poles 

Res F[l^-^(9 + iir,9 ,9 1 ...,9 k ) = i F^"** (9 U . . . A), (2.11) 



k+2 

9o = 0o 



Res F£f^*(0 o + iMo,0i...A) = -iJlS^^FP^^, ... ,9 k ). (2.12) 



?0 = »0 



1=1 



Here ft = (a, c) is the anti-particle of fi. The equations (|2.1ip and (|2.12p are in fact not 
independent from each other, as solutions to (|2.12p may be obtained from the solutions to 
(|2.1ip by employing the first two form factor equations. This is the reason why later on we will 
concentrate on solving (|2,lip . 



2.1 Minimal form factors 

In order to solve equations (|2.8p - (|2.12p it is common practise to adopt a recursive approach 
whereby one starts by computing the two-particle form factors (the zero-particle and one-particle 
form factors are constant for spinless fields, with Fj~ = (0|T|0) := (T)) and then uses equations 
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(|2.1ip - (|2.12p to obtain higher particle form factors iteratively. A minimal solution (e.g. with no 
poles on the physical sheet) can be obtained by solving the two particle versions of equations 
(123]) and $2JB, that is 



f2^ 2 (9) = S^ 2 {e)Fl^{-e) = F^T 1 (2*™ - 9), V Ml , M2 (2.13) 

which, combined with the requirement that the two particle form factor must have simple poles 
at 6 = in and 6 = ii:(2n — 1) and that the residue at the first pole must be 

Res F^ m {6 + m) = iF^:=i{T), (2.14) 
6 = 

gives the general solution [24j 

f T\ih^iq\ = (7~) sin(|) i? m"in lfe ( 6> ) (2 15) 

2n sinh ( M2( Cl -c2)-i) + ^ ginh (HtMdH) F ^ (<7r ) ' 

where ci 5 2 represent the copy number associated to each of the particles as defined in ()2.7p . 

3 Introducing the models and the A-sum rule 

In |24j a general expression for the two-particle form factors of the twist field was found that 
applies to all integrable QFTs. This expression was then specialised to the Ising and sinh-Gordon 
models and checked against the A-sum rule [30] . In its original form, the A-sum rule can be 
expressed as 

1 f°° 

A° = -—— drr(O(r)O{0)) c , (3.16) 



2(0} Jo 

where the subindex c in the two-point function stands for "connected" , meaning that the product 
of the vacuum expectation values of the fields has been subtracted. is the trace of the 
energy momentum tensor of the model. Eq. f|3. 16j) provides an expression for the conformal 
dimension of a primary field. The remarkable fact is that the dimension is given in terms of a 
two point function involving the field O which represents the counterpart of the primary field 
in the perturbed (massive) model. Therefore (|3.16|) allows us to extract information about the 
underlying CFT from the two-point function of a massive theory. 
A slightly more general version of (|3.16p was employed in [37] 



i r°° 

A ° {r ° ) = ~ 2(0)J r drr ^°^^ ( 3 - 17 ) 

so that taking ro = we recover (|3.16p whereas for larger values of ro we are now able to trace 
changes in the value of A^(ro) along the renormalisation group (RG) flow, that is as we move 
from low energies or ro large to high energies or ro = 0. Observing such intermediate behaviour 
is particularly interesting for models where the RG-flows approach the vicinity of more than one 
critical point, as the ones we will consider below. 

Taking O = T, employing (|2.5p and performing the integration in r (|3.17p becomes 



d6i . . . d6 k (1 + r E) e 
2 (T) ^ ^ J " J k\(2ir) k 2E 2 



-rnE 



-oo 



xF^-^(9 1 ,...,9 k ) (fJ Im *(0 1 ,...A)) , (3-18) 
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where E stands for the sum of the on-shell energies E = Y^,i=i m M- cos (^«) [3Z]- I n those models 
for which the form factors of G and T are known, one can in principle identify the conformal 
dimension A' by computing (|3.18p and taking r$ = 0. This dimension is known a priori to be 

Ar (0 ):=A'-=± („-!), (3.19) 

so that evaluating (|3.16p effectively allows for a consistency check of the form factors of T ■ Such 
a check was carried out successfully in [23] for the Ising and sinh-Gordon models. There only 
the two-particle contribution to the expansion of the two-point function was considered and very 
good agreement with the predicted value was reached (the agreement was exact for the Ising 
model). The reason why the two-particle approximation works so well is that the expansion 
(|2.5|) is in fact rapidly convergent as a function of the particle numbers involved in the form 
factors. Therefore, in general the two-particle contribution is the more major and for many 
models other contributions are almost negligible. 

Although this is a very useful feature which has been widely exploited for integrable models, 
it is inconvenient if one wants to use (|3.18p to test higher particle form factor solutions. There 
is however a way to overcome this obstacle, that is, by considering models for which the two 
particle contribution is far from providing a good picture of the ultraviolet behaviour of the 
theory. In the upcoming subsections we will describe two models which have precisely this 
feature: the roaming trajectories model and the SU (3)2-homogeneous sine-Gordon model. 



3.1 The roaming trajectories model 

The first theory we want to investigate here is the roaming trajectories (RT) model [38] . This 
is a model with a single particle spectrum and no bound states which is closely related to the 
sinh-Gordon model. The model is characterized by the two-particle S-matrix, 

S(0)=tanh- (0-0 o - - J tanh - (0 + O - y J , GR. (3.20) 

On the other hand, the sinh-Gordon S-matrix \39\ BO] is given by 

, v tanh i (0 - ^) 

*<»>==4t»T# <3 - 21) 

It is easy to see that the S-matrix (|3.20p can be obtained from (|3.2ip by the replacement 

B^l-^. (3.22) 

This relationship implies in particular that computing the form factors of the sinh-Gordon model 
and setting B to the value (|3.22|) gives the form factors of the RT-model. 

The roaming trajectories the model's name refers to emerged in the computation of the 
effective central charge c e g (r) within the thermodynamic Bethe ansatz approach [41|, EE2"] carried 
out in [38]. For massive QFTs it is expected that the function c e ff(r) "flows" from the value 
zero in the infrared (large r) to a finite value in the ultraviolet (small r). For many theories, 
including the sinh-Gordon model, the constant value reached as r — > is the central charge 
of the underlying conformal field theory associated to the model. In this case, that theory is 
the free massless boson, a conformal field theory with central charge c = 1. Therefore, in the 
sinh-Gordon model, the function c e g(r) flows from the value zero to the value 1 as r decreases. 
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Crucially, when the same function c e s(r) is computed for the RT-model it shows a very 
different behaviour. It still flows from the value to the value 1, but it does so by "visiting" 
infinitely many intermediate values of c giving rise to a staircase (or roaming) pattern. The 
values of c that are visited correspond exactly to the central charges of the unitary minimal 
models of conformal field theory 

Cp = l- . , 1V with p = 3,4,5... (3.23) 
p(p + l) 

Another observation made in [38] is that the size of the intermediate plateaux that the function 
c eff( r ) develops at the values (|3.23p is determined by the value of 9q. For 6q = there is a single 
plateaux at c = 1, thus the usual sinh-Gordon behaviour is recovered, whereas the plateaux at 
(|3.23p become more prominent as 9q is increased. In the limit 9$ — > oo a single plateaux at 
c = \ remains which reflects the fact that the 5-matrix (|3.20p becomes -1 in this limit, hence 
the model reduces to the Ising field theory. This interesting limit behaviour was studied in [33] 
within the form factor approach. 



3.2 The S'?7(3)2-Homogeneous sine-Gordon model 

The second model we want to study is the 5 , C/(3)2-Homogeneous sine-Gordon (HSG) model. 
The model is just one of the simplest representatives of a large class of theories first named in 
[33] whose spectrum [33J [35l [36] , S-matrix [37] , form factors [3HJ [39l [37J [50] and thermodynamic 
properties [5U [52l [53] have been extensively investigated over the last two decades. The HSG- 
models are very interesting theories, as they include a number of distinct features rarely found 
for integrable models: they posses both unstable particles and bound states in their spectrum 
and their S'-matrices are generally non-parity invariant, that is S a b(0) ^ Sb a (9) for a ^ b. In 
particular, the S'C/(3)2-HSG model contains two particles, which we will label as + and — . They 
are self-conjugated and interact with each other by means of the following S-matrix 

S±±(9) = -l, and S± T (9) = itanh^ U±a- % ^\ . (3.24) 

Thus particles of the same species interact with each other as free fermions, whereas particles of 
different species interact by means of parity-breaking S-matrix which depends on a free param- 
eter a. These S-matrix amplitudes have a pole in the unphysical sheet (that is ^(9) S (— 7r,0)), 
with real part given by ±cr. Such type of poles are a signature of the presence of unstable 
particles in the spectrum. 

The scattering picture is that particles + and — interact with each other by creating an 
unstable particle, whose mass and decay width depend on the parameter a through Breit- 
Wigner's formula [53]. More precisely, for \a\ large, the mass of the unstable particle can be 
approximated by me' 17 '/ 2 , where m is the mass of the stable particles [37] . Therefore, the limit 
a —> oo corresponds to an infinitely massive unstable particle, that is a particle that can not be 
formed at any finite energy scales. At the level of the S-matrix we find that lim -_ 5 . 00 S± t ^(9) = 1, 
that is, the model reduces to two non-interacting copies of the Ising field theory. This property 
is very useful as a consistency check in form factor calculations. It implies that when a — > oo the 
form factors of any field should reduce to those of the Ising model, which are generally known. 

As for the RT-model described before, the effective central charge of the St/(3)2-HSG model 
also exhibits a staircase pattern, albeit with only two steps (at most) [51] . The same structure 
was found for Zamolodchikov's c-function and the conformal dimensions of certain local fields 
[55[ [37] . In this case the appearance of steps is directly related to the presence of the unstable 
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particle and its mass. There is only one step if a = in which case the unstable particle's mass 
is of the same order as that of the stable particles and a second step emerges if a 7^ whose 
onset and length are related to the precise value of a. All these features have been analysed in 
detail in |51|, 137] . In section 6 we will see that the conformal dimension of the twist field (|3.18|) 
is no exception to this behaviour. 



4 Twist field form factors for the RT-model 

We will start our analysis by considering the simplest of the two models described above, in 
terms of its particle spectrum. The two-particle minimal form factor of the sinh-Gordon model 



f2L U (*) = eip 



, 00 (ftsinh^sinh^^ 

-2/ — ±— -J_cosht(n + - 

Iq t smh(nt ) cosh | 



t sinh(nt) cosh | \ ir 



(4.25) 



was first obtained in [24] and can be easily rewritten as an infinite product of ratios of Gamma 
functions. The explicit expression can be also found in [24] . 
It is natural to make the following ansatz, 

k ^ r l n (B_) 

Fl( Xl ,...,x k ) = H k Q k (x u ...,x k )ll — mm Xj _ -r, (4.26) 

. yXj, OLX j J yX j QtOC<i J 

where we have introduced the new variables x% = e n and a = e n so that, for example 



A similar ansatz was already used in [56J in a different context. We use the simplified notation 
FT(xi, . . . ,x k ) to represent the A;-particle form factor of particles all of which live in the same 
copy of the model. The functions Qk{x\, Xk) are symmetric in all variables and have no poles 
on the physical sheet. H k are rapidity independent. 

The ansatz (|4.26p is reminiscent of the solution procedure that is traditionally used in the 
original form factor programme (see e.g. [S7] where the sinh-Gordon model was studied). This 
ansatz is useful as it isolates the pole structure of the form factors in the product. Provided that 
Q k (xi, ...,x k ) are analytic functions, symmetric in all variables, then it automatically satisfies 
equations (|2.8p and (|2.9|) . For k = 0, the condition (|2.14p implies the normalization Hq = (T) 
and Qq = 1. 

Once the ansatz (|4.26p has been made it remains to identify the functions Q k (xi, . . . ,x k ) 
and the constants H k . In the sinh-Gordon model symmetry considerations imply that only even 
particle form factors are non-vanishing, so that our first new results would correspond to the 
k = 4 case and k will always be an even number. We therefore turn to solving equation (|2.1ip . 
which we can now rewrite as 

lim (0 O - o )F k+2 (ax o ,x Q ,x 1 , . . . ,x k ) = iF k (xi, . . . ,x k ), (4.28) 
o ->0o 

where xq = e"». 

In order to turn the equation (|4.28|) into an equation for the functions Q k (xi, . . . ,x k ) and 
the constants H k the following identity will be needed, 

T |n qx T |n x _ (x -Xi)(ax - Xj) 
Xi Xi (a/3 x Q - Xi){/3x Q - Xi) 
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where f3 = e 2™ and B is the coupling constant that appears in the sinh-Gordon S'-matrix (|3.2ip . 
This identity can be easily derived from the Gamma function representation of the minimal form 
factor [23| • 

Substituting the ansatz (|5.52p into (|4.28|) and simplifying we obtain 
2 sin - a k+2 

H k+2 = H k and Q k+2 (ax ,x ,xi, . . . ,x k ) = xlP k Q k (xi, . . . ,x k ). (4.30) 

nF mL ("0 

where 

k 

Pk = {x a ~ a 2 x )(x - ax b )(x c - af3~ l x q )(I3x - x d ) (4.31) 



a,b,c,d=l 

k 

{-a) k £ (-a 2 xo) fc - a (-a" 1 xo) fc ^(-ar 1 ^o) fc " c (-^o) fe -^)^ fe) a( fc ) C Ti fe) , 

a,b,c,d=0 



(k) 

and a \ is the i-th elementary symmetric polynomial on k variables x\, . . . ,x k , which can be 
defined by means of the generating function, 

k k 

Y,x k - l a\ k) =H(x l + x). (4.32) 

i=0 t=l 

The equation for H k can be easily solved to 

k 

( 2sin-a 2 \ 2 kf k ^ ,_ v 
Hk =[ w T\n r J "^Hn (4-33) 

whereas equations for the polynomials Q k (xi, . . . ,x k ) will need to be solved on a case by case 
basis. Unfortunately the solutions get very involved very quickly. There are three main reasons 
for this: 

• The degree of the polynomial in the denominator of (|4.26p is much higher than would 
be the case in the standard form factor programme. Since the twist field is spinless, the 
degree of such polynomial must equal the degree of the polynomial Q k (x\, . . . , x k ) and this 
means that its degree will be very high for relatively small values of k. As an example, for 
the RT-model we will see later that the degree of Q2(x\,X2) is just 2, but the degrees of 
Qi{x\, X2, X3, X4) and Qe(xi, . . . ,xq) are 12 and 30 respectively. 

• The polynomial P k is a very complicated function in terms of elementary symmetric poly- 
nomials, which again complicates the solution procedure and makes it very difficult to 
identify any patterns as k is increased. 

(k) 

• The reduction properties of the elementary symmetric polynomials cr- are much more 
involved for the twist field than in the usual form factor programme. In general, 

af +2 > = of + (1 + ajzooffi + (4-34) 

where at +2 ^ is an elementary symmetric polynomial on the variables axo, xq, x\, . . . ,x k 
and a t^ , cr t—\ 1 ® 1— 2 are elementary symmetric polynomials in the variables x\, . . . ,x k - We 

(k) (k) 

will also adopt the conventions a\ = for i < and a = 1. The usual reduction 
properties are recovered for n = 1 or a = -1. 
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The polynomial Q 2 (xi,x 2 ) can be easily obtained by setting k = in (|4.30p which gives the 
equation, 

Q 2 (ax ,x ) = xl, (4.35) 

(21 

There are actually two combinations of elementary symmetric polynomials of two variables a\ 

(2) 

and a 2 that correctly reduce to the identity above. The most general solution is 

Q2{Xi,x 2 ) = a a\' + ^2-^2(^1,^2), (4.36) 
with f^2 an arbitrary constant and 

/ (2) \ 2 

K 2 (x 1 ,x 2 ) = a~ l af - I ^— J , (4.37) 

the kernel of equation (|4.35p . that is the most general order 2 polynomial on the variables xi,x 2 
which solves 

Q k+2 (ax ,x ,xi, . . . ,x k ) = 0, (4.38) 

with k = 0. Substituting (|4.36p together with H 2 in (|4.26|> it is easy to see that (j2.15[) is only 
recovered for p l = p 2 = c± = c 2 = 1 if we choose Q 2 = 0. Hence we have fixed the constant 
above and can now go on to compute the four particle form factor. 

Solving now for Qi(x±, x 2 , X3, X4) we find that the most general solution to (|4.3U|) takes the 
form 

Q4,(xi,x 2 ,x 3 ,X4) = 04 [a 2 + 702(03 + 0-10-4) + Serous + r\a\a\ + io\o± 

+A0-10-30-4 + paj\ + tt^K^xi, x 2 , X3, X4), (4.39) 

where we have abbreviated = a{. The constants 7, 5, 77, £, A and p are fixed functions of n, 
whose explicit form is given in appendix A. The function K&(xi, x 2 , X3, X4) is the most general 
order 12 polynomial on the variables xi,x 2 ,x 3 and X4 that solves the equation (|4.38p and Q4 is 
an arbitrary constant. The function K^{x\, x 2 , x 3 , £4) has the form 



K±(xi,x 2 ,x 3 ,x 4: ) = Aa\ala\ + B(a\al + a\a\ + a\ala^) + Ca 1 a 2 a 3 (a\ + a\i 

A A O A /O 

+ Da 2 a<± + 0-i0- 4 + 0-3 + -E0-10-20-30-4 + Fa^a^a^ + Ga 2 a^[a 3 
+ Ha 2 a\ + /010-3O4 + J0-4, 



3 + o?o 4 ) 

+ 010-4) 
(4.40) 



where the constants are given in appendix A. 

We have also computed the most general polynomial Qq(x\, . . . , xq) which solves ()4.30p with 
k = 4. The solution is an order 30 polynomial on the variables x±,...,xq and too cumbersome 
to be reported here. For O4 = (we will see below why this choice is sensible), Q§{x\, . . . ,xq) 
depends once more on a free parameter Qq, which as above acts as coefficient to the function 
Kq(x\, . . . , xq) which satisfies the same equation (|4.38p above. 

Therefore, a structure seems to emerge where the most general 2£>particle form factor de- 
pends on k free parameters. A similar structure was found when studying the boundary form 
factors of specific fields in the ^-affine Toda field theory [251 ES], although no physical inter- 
pretation for the result was provided there. A more thorough analysis of solutions to equations 
of the form (|4.38|) was carried out in [BO] for the case a = — 1 and the field TT. 
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4.1 Identifying the twist field form factors 

In this section we would like to argue that choosing O4 = in (|4.39[) corresponds to the specific 
twist field we are interested in. The general solution (|4.39j) is a one-parameter family of solutions 
characterized by the choice of the constant ^4. Given the usual assumption that the space of 
fields in a local QFT is linear, we expect that the form factor of a linear combination of fields is 
a linear combination of form factors, that is, in general 

F fc 0l+n0a (xi, . . . ,x k ) = F&ixu . . . ,x k ) + Oif 2 (x!, . . . ,x k ) (4.41) 

and therefore the solution (|4.39p must describe the form factors of a linear combination of local 
fields (as would the solution (|4,36p ). Since we are interested only in one very particular field, 
the twist field 7~, we must find a suitable mechanism that allows us to select the particular value 
of corresponding to the four-particle form factor of the twist field. 

An interesting way of identifying the form factors of the twist field is to use the form factor 
cluster decomposition property, which has been studied for various models in the past [ST1 ESI 
163"! I49j and analysed from a more general point of view in [30]. It is a factorization property of 
form factors which, for the four particle case, can be expressed as 

lim fZ(kxukx2,x 5 ,x a ) oc F^(x uX2 )F?(x a ,x A ). (4.42) 

In general, the fields 71 and T2 on the r.h.s. may not necessarily correspond to the same field as 
the form factor on the l.h.s. A notable example of this is the model studied in [49] and the form 
factors of the field TT studied in |60| . In [3U] it was argued that for theories without internal 
symmetries, the cluster decomposition would be a consequence of the decoupling of right- and 
left-moving modes in the conformal limit and would hold for any field whose counterpart in the 
underlying conformal field theory is a primary field. 

Given that the twist field does certainly correspond to a primary field in the underlying 
conformal field theory we expect a factorization of the type (|4.42p . Imposing (|4.42p is in fact 
sufficient to select a single value of in (|4.39p . Indeed, if we carry out the cluster limit in 
(|4,42p for the general expression (|4.39p and we call <Ji = (Ji{x\,x 2 ) and &i = ai{x^,x^) we find 
that 

lim FAkxi,kx2iXz : X4) ~ |0"2<5"2 + £lAAa\a\ + B(a 2 a\ + crf<7 2 ) + Da2a 2 )] 

x mm \x 2 > mm \x A > ,^ 

(x\ — ax2){x2 — axi)(xs — ax\) (24 — ax^) 

Clearly, this expression factorises if and only if O4 = 0. In that case, we recover exactly (|4.42l) 
with 71 = T2 = T ■ We will therefore choose Q4 = as our twist field solution. 

If we had chosen to use the cluster decomposition property to fix the constant m (|4.36p 
we would have found 

lim F 2 r («;xi,X2) oc fi 2 , (4.44) 

K — ^OO 

so that our choice ^2 = guarantees that lim K _ i . 00 F^kxi, x 2 ) oc FJfJ = 0. 

In general, it appears from our two-, four- and six-particle form factor solutions that for every 
a G Z + there exists a field /C 2 a whose form factors solve (|2.8p - (|2,12p and have the interesting 
property that 

if 2 " = F^^{B u e 2 ) = ■■■ = F^-^-^Ox, . . . , 6 2a _ 2 ) = 0, (4.45) 
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consequently i 7 ^^" . . . , O20) solves ()4.38p for k = 2a, that is, it has no kinematic poles. 

At this stage we can only speculate about the nature of the fields K,i a - From the a = 2 
example, the cluster decomposition property suggests that the field /C4 does not correspond to 
a primary at conformal level, since at least its four particle form factor does not factorize under 
clustering. Furthermore, given that the form factors of all fields YZia axe solutions to the twist 
field form factor equations they must be twist fields of some kind. Finally, their non-vanishing 
form factors involve even particle numbers, which points to a particular kind of symmetry. One 
may think of linear combinations of composite fields such as OT or tp 2 7~ etc. As a future project 
it would be very interesting to identify the precise nature of the fields K,2a- 



5 Twist field form factors for the S£/(3)2-HSG model 

We turn now to the second model whose twist field form factors we wish to investigate. The 
complexity of the model is increased by a number of features, notably the fact that its spectrum 
has two particles and the presence of the free parameter a in the S'-matrix (|3.24p . The solutions 
to the equations (|2.13|) with the S-matrices (|3.24p are 

^ (*) = "isinh (^), (5.46) 



and 



= A(n)e±t+^ exp [ f X \ + ' I ■ (^T) 

t sinh(nt)cosh(t/2) 1 ' 



-00 




with A(n) given by the limit, 



2 P +2+n in 2n(4p + 3 + 2n\ sn /4p + 5 + 2n\ 8n -i-r 1 ( 4 „ ) 
A(n) = lime ^ 4 . — — i 22 — ' 

4 " (5.48) 

The solution (j5.46p is nothing but the Ising model solution first obtained in |24] . as we would 
expect from the first S-matrix in (|3.24D . 

The form factors (|5.47p can also be expressed in terms of an infinite product of Gamma 
functions 

00 r ( ik+3+2n\ 2 t / -2w+ik+l+2n\ r /2w+4k+l+2n\ 

= A(n)e^+^ TT { ^^h ( 4n ] ( An K (5-49) 

mmV ' y ' -I- -I- r j Ak+l+2n \ 2 -r I -2w+4k+3+2n \ r ( 2w+ lfc+3+2n \ ' v ' 

k=0 L { in ) 1 I 4^ ) L \ 4^ ) 

with w = n + i{9 ± cr)/ir. 

The function A(n) defined above would seem a strange choice of normalization. The moti- 
vation for it is to ensure that the following minimal form factor relations 

I 9 j in 

*&V)*&(0 + -) = ± J"" y (5.50) 

blml \ 2n "1" An) 

hold, without involving complicated constants. In particular, A(l) = e~ G / 7r e~ i7r//4 2 1//4 , where G 
is the Catalan constant that appears in the normalization of the form factors of the one-copy 
model [64^ 148] . It is worth noticing however that with respect to the latter normalization our 
minimal form factor at n = 1 is multiplied by the extra factor e _l7r//4 . 
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Once the two-particle form factor and minimal form factor have been computed the basic 
monodromy and pole structure features of the form factors are fixed so that higher particle 
form factors can be constructed in terms of the solutions already found. Let us introduce the 
following notation: 

£ m 

F i+m ({x}j; ) := ' ' ' + ~ ' ' ' ~ (x\, . . . , x £ , x e+1 . . . x e+m ), (5.51) 

This represents the I + m-particle form factor of the twist field with t particles of type + and m 
particles of type — living in one particular copy of the model. For the model under consideration, 
we will make the following ansatz 

pT\++ /x±\ 

F i+m ({x}f;{x}-) = H+-Qj-J{x}j;{x}-) ]J 



\ Xj GiX n ) \ X n CXXj } 

i<i<j<e y 3A ■> % > 



£ l+m T?T\ ix^\ 

n n F ^~i-) n , - ( 5 - 52 ) 

X n \ Xj GiX j ) \ X j GiX 'ij ) 

i=Xj=i+l ■> e+l<i<j<£+m y J,K ■> ' 



In terms of the new variables X{ we can rewrite for example 



F^Qi-Oj) = (5-53) 

Xj 

It is easy to check that, the ansatz (|5.52|) automatically satisfies equations (|2.8p and (|2.9p 
provided that the functions Q~^+ m {{x}~£~ '; ) are separately symmetric in both sets of variables 
and have no poles on the physical sheet and Hf~ are rapidity independent. Notice that there 
are kinematic poles associated to pairs of + and — particles, but not to the combination +— , 
as the two particles in the model are self-conjugated (their own antiparticle) . The ansatz (|5.52l) 
is reminiscent of the solution procedure used in |48t H9] where the form factors of local fields of 
the present model were also studied. 

Once the ansatz (|5.52p has been made it remains to identify the functions Q^ m ({^}/; {x}^) 
and the constants Ht~- A useful benchmark that can be employed for this model is the fact 
that whenever m = or £ = 0, the resulting form factor must be the ^-particle or ?n-particle 
form factor of the Ising model, respectively. This relationship with the Ising model, combined 
with the kinematic residue equation (|2.1ip also implies that only form factors with both I and 
m even will be non- vanishing. 

Substituting the ansatz (|5.52|) into (|2.1ip we obtain the following recursive relations for 
Qt+ m ({x}t; {x}m) an d tne constants H^~, 

3£ — m i o am rj/j _ i i 

+ _ _ a^- +2 e-^r2^- m+1 sin2 + _ 
H e+2,m = ,-,71++/ s ~ ^<,m' (5.54) 

nF mL 

and 

Qt+2+m( ax 0> x 0,{ x }t'>{ x }m) = P £m( x od x }t'd x }m)Qt+m({ x }t'd x }m), ( 5 - 55 ) 

with 

£ „_ . m 

^^,{^;MJ=«T-V Y, "\-a 2 ^V+ -+£(-V5e£xo) m -*^, 

i,j=0 k=0 



(5.56) 
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where , oT are elementary symmetric polynomials on the variables {x}/ and {x}~, respec- 
tively. To simplify notation, in (|5.56p and (|5.59p we have dropped the explicit variable depen- 
dence of the symmetric polynomials. 

In the ansatz (|5.52|) we have chosen a particular ordering of the particles with type + 
appearing first and type - last. Of course this ordering can be changed by employing the 
first form factor equation (|2.8p . Alternatively, we could have worked with the form factor 
Fe+ m ({x}j ; {x}+) where we now have t particles of type - first, followed by m particles of type 
+. For this ordering, the recurrence equations above become instead 

a f+2 e -^2 2 ^ m+1 S in^ , , s 

H I+2,m = "TV— n ~ Hi+, (5.57) 

nF min (**") 



Q,X+m(oixo,xo,{x} e ;{x}+) = P+(x ,{x}, ; {x}+ )Q , + m {{x}„ ;{x}+), (5.58) 



and 
with 

t 



m— fc _+ 



m ij=0 fc=0 



(5.59) 



From the definition (|5.52p and equations (|2.8p and (|2. 13|) it is easy to show that 

BfaQfcii*}* tom) = H-+Q~+({x} m ; {x} t ), (5.60) 
which provides a useful relationship between the solutions of (|5.58p and those of f[5.55]> . 

5.1 Solutions to the recursive equations 

Given the structure of the ^-matrix (|3.24p we know that form factors involving only particles 
of type + or only particles of type - should equal the form factors of the Ising model. We will 
therefore split our solutions into Ising model solutions and solutions involving particles of both 
types. 

5.1.1 Ising model solutions 

For m = in (|5.54p - (|5.55p or equivalently m = in (|5.57p - (|5.58p the equations reduce to the 
form factor equations of the Ising model. That is, 

a f+2 2 2*+l sin2 r 

and 

Qe +2 (ax , x , {x}i) = Pi(x , {x}e)Qi({x}e), (5.62) 

with 



I 

P e (x ,{x} e ) = a e x i + 2 a e Y, ("^) (-c^xoY^aiaj. (5.63) 

i,j=o a 
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Interestingly even for the Ising model, these equations are not easy to solve and the solutions 
for Qe({x}e) become very cumbersome beyond £ = 4. The first few solutions are, 

Q2(xi,x 2 ) = a~V 2 , (5.64) 

n ( \ -12(2 Pi(«)o*io-3 . Pi(a)(l + a 2 ) 2 a 4 \ , > 

(^4(^1,^2, £3,^4) = a °4 °2 1 q h (5.65) 

\ a a 6 ) 

3 / 2 2, Pl(a0 2 O-l<7i<7 5 p 2 (a)0-l0-20-4<7 5 

\ or a 

(l + a 2 ) 2 p 3 (a)<rfc§ _ pi(a)g 3 (o-i^ + ^^s) 
a 4 a 
, P2(a)pi(a)V|a 6 (l + a 2 f Pl {a) (a\ + a|fj 6 ) 



a 5 a 3 



P 2 (q)pi(q) P4(Q)P3(«)q"iQ-5C r 6 p 2 («) Pi(a) P3(«)^| 

a 7 a 9 
Pi(«) 2 P3(a)o"3 (040-5 + criO- 2 cr 6 ) P2(a)pi(a) 2 p 5 (a)o-20-4CT6 



a 4 



Pi{a)p3{a) 2 (0-20-5 + afa^e) 

a 5 



(5.66) 



with 



Pi 1 




9 

= 1 + a + a , 


p 2 l 




= 1 — a + a 2 


P3l 




= 1 + a + a 2 + a 3 + a 4 , 


P4l 




= 1 — a + 3a 2 — a 3 + a 4 , 


P5l 




= 3 + 2a + 4a 2 + 2a 3 + 3a 4 



(5.67) 

However the twist field form factors are already know for the Ising model. They were computed 
in [28] not by solving equations (|2.8|) - (|2.12p but by using the special free fermion features of the 
model. It was found that 

Fl(0 1 ,...,e t ) = Pf(KW), (5.68) 

where Pf(K^) = \J det(K^) and is an £ x £ matrix whose entries are given by = 
-^2 )• Comparing to our original ansatz we have the remarkable identity 

i 

Qt{{x}i) = H^V^K 1 ) Y[(xi - a Xj ){ Xj - axt). (5.69) 

i<j 

Bringing the r.h.s. of (|5,69p into the form of a combination of symmetric polynomials is highly 
non-trivial for £ > 4. In particular, for I = 6 it yields the result ()5.66p . It would be nice to 
develop a general proof of (|5.69|) . 
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5.2 Solutions involving both particle types 

Starting with the two particle solutions (|5.64p we find the following new four particle form factor 
solutions 



Q^ +2 (xi,x 2 ;z3,£4) 
Q^+ 2 (xi,x 2 ;x3,x 4 ) 



a (J, a 



a <Tn a 



l + a 
fa 



o{ a 1 + a 



2 ' 



l + a 



2 > 



(5.70) 
(5.71) 



where are symmetric polynomials on the variables {xe™}^ m . 

Going beyond four particles is rather difficult, but because of their relationship to the form 
factors of the Ising model, it is possible to find closed formulae for certain types of form factors. 
For example, when 1 = 2 and m is general. In this particular case the form factor equations 
become simply 



Qi+m( ax o,x ;{x} m ) = x\ m ^(-V / ae»3;o) m k a k Q m ({x}] 



k=0 



and 



Q2+m( ax o,x ;{x}; 



„2 m 



Y,(-^x ) m - k a+Q m ({x}; 



(5.72) 



(5.73) 



fc=0 



where Q m ({x}m) is the Ising model solution given by (|5.69|) . Particular solutions to (|5.72|) and 
(I5.73P take the form, 



in 

(xi,x 2 ;{x} 

i 



a 2 



2k 



m — 2 

m— 1 I 2 

a 2 aT x-^ a 2k+i 



+ \k 



(-2 ) 



-E 



Qm.({^}t 



(5.74) 



and 



'2+m 



(xi,x 2 ;{x}; 



a: 



-^E 



(7 



2A- 



m-2 

i ^ °"2A:+1 



fcS (^2 f 



Qm(W+). (5.75) 



They provide closed solutions to the equations (|5.T2j) - (|5.T3j) valid for any values of m. Unfortu- 
nately, this is not enough to conclude they are fully consistent with all form factor equations. 
What we mean is that the relation (j5.60j) must also hold, which means that for example the 
solution Q^+^xi, x 2 ; {x]\ ) constructed above, must also solve the form factor equation satis- 
fied by QXZi{\_x\\\x\^x-i) (up to constants). This imposes a set of further constraints on the 
solutions to (|5.72|) and (|5.73|) . Let us consider an example. 
From (|5.74p we find 



Q 2 + 4 (xi,x 2 ;{x} 4 ) = a 



(a 2 -) 2 + a~ 2 at + a+ - («r 2 a+ + a+) 



Q^{x}+). (5.76) 



This function solves (|5.72p . however it does not solve the equation for Q~l+2{{^}t ; x±, x 2 ) which 
can be obtained from (|5.59p with £ = 2, m = 4. If we solve that equation, we obtain a completely 
different solution. Therefore (|5.76p is not a consistent solution to all form factor equations. As 
we studied in detail for the RT-model, we can generally add an extra function to any solution, 
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as long as that function is in the kernel of the equation we are trying to solve. In our case, this 
means that we can always add to ()5.76p any function K 2 ^ A (x\, x 2 \ {x}^) which satisfies, 

K 2+4 ( ax o, x ; {x}j) = 0. (5.77) 
The most general solution to this equation, up to a multiplicative constant is, 

cr(l + a) z 

x (a 3 afa%(Tt -a(l + a + a 2 )((a^) 2 + (affa^) + (1 + a) 4 (l + a 2 )a^at)) . (5.78) 
Solving for ± 2 ({x}j[ ; xi, £2) we find that 



(<r 2 ) 2 + <r 2 a$ + a\ - {a 2 af + a+) 



Q 4 ({x}+) 



+K 2 + 4 (x 1 ,x 2 ;{x}+), (5.79) 

and from equation ()5.60p it follows that 

Q 2 " 4 f 4 (xi,x 2 ;{x}J) = a~ 2 Qj~ 2 ({x}j;x 1 ,x 2 ). (5.80) 

Therefore, in general, the solutions (|5.74|) and (|5.75p need to be modified by adding some function 
in the kernel of (|5.72p or (|5.73|) which is consistent with (|5.60p . 

6 Numerical results 

In this section we want to provide numerical evidence that the twist field form factors computed 
thus far do indeed correspond to the correct twist field. Our method is to check numerically the 
form factors of the two models presented in sections 13.11 and 13.21 against the A-sum rule. For 
both models, we considered only the two and four particle terms in the expansion (|3,18p . 

For each model, we have adopted a different numerical approach: For the RT-model we have 
only evaluated (|3,17p directly in the ultraviolet limit r$ = 0. This is because the numerical 
recipe used was particularly slow. We employed a truncated version of the infinite product of 
Gamma functions given in [23] to evaluate the minimal form factor (|4.25[> . We also employed a 
Newton-Cotes method and a Montecarlo simulation which uses the Vegas algorithm to perform 
the integrals in the two-particle and four-particle case, respectively. 

For the HSG-model though we have been able to consider a wide range of values of tq by em- 
ploying a very precise, piece- wise polynomial interpolation of the functions F^J^ (9) (6)* 
which dramatically reduces the running time of our programme. In this case we have carried 
out the integrals by means of the Vegas algorithm. 

The numerical results are shown below. 



6.1 Numerical results for the roaming trajectories model 

As explained in section I3TT1 the function c c fr(r) or effective central charge of this model exhibits 
an infinite set of plateaux between r = and r — > 00. A similar type of behaviour is expected 
for A^(ro) as ro is varied. Here we have only considered ro = 0, however our results still allow 
us to identify two values of A 7 ", that is the value obtained in the two-particle approximation and 
the value obtained in the four-particle approximation. Each of these values agrees with what 
we would have expected for the first two plateaux of the function A 7 "(ro). 
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The precise location of the plateaux can be easily predicted by combining f|3.23[) with (|3. 19j) . 
The first two values of the central charge correspond to C3 = 1/2 and C4 = 7/10, respectively. 
Inserting these values in (j3. 19|) we obtain a value of A 7 " for each central charge and each value 
of n. 

The two particle contribution takes exactly the same form as for the sinh-Gordon model and 
was given in [23]. Evaluating it for 6q = 20 we obtain the values listed below 



n 




Al(0) 


2 


0.03125 


0.0312548 


3 


0.0555556 


0.055676 


4 


0.078125 


0.0785953 


5 


0.1 


0.101033 


6 


0.121528 


0.123257 


7 


0.142857 


0.145351 


8 


0.164062 


0.167351 


9 


0.185185 


0.189277 


10 


0.20625 


0.211143 



Table 1: Two particle contribution to the conformal dimension in the RT- model. The second 
column shows the exact values of the conformal dimension of the twist field corresponding to 
central charge C3 = 1/2. The third column shows the numerical values of the same quantity in 
the two-particle approximation for #0 = 20. 



Employing the four-particle form factors of the energy-momentum tensor obtained in [57] and 
and our solution (|4.39|) with = 0, the four-particle contribution is given by 

T, n , sin(^)cosh(go)|F^ in (ivr)-FQ n (m)| 2 



Ai (0) 



1536vr 3 



U X a 2 (TzQi{x 1 ,X2,X^X i ) n i^in^^mmV^j; 

de 1 de 2 de 3 de 4 — -. ^ (6.8i) 

4 



00 



yicos(£) -cosh^J j (Ecosh(W 

where <7i,a 2 and C3 above represent elementary symmetric polynomials in the variables e di 
with i = 1,2,3,4. The values of (|6.8ip for different values of n are given in Table 2. Both 
Table 1 and 2 show relatively good agreement between the values predicted by the theory and 
those numerically obtained. The difference between the theoretical and numerical values is 
considerable for some of the results in Table 2, specially as n is increased. However it is always 
within the standard deviation of the computation. 

6.2 A-sum rule for the <W(3) 2 -HSG model 

From (|2.15p and the two particle form factor of the energy-momentum tensor for the thermally 
perturbed Ising model 

F 2 6|±± (0) = -2mm 2 sinh(^), (6.82) 
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n 




ARO) 


2 


0.012500 


0.013086 


3 


0.022200 


0.022169 


4 


0.031250 


0.028611 


5 


0.040000 


0.042555 


6 


0.048611 


0.047566 


7 


0.057143 


0.057996 


8 


0.065625 


0.064736 


9 


0.074074 


0.072281 


10 


0.082500 


0.068762 



Table 2: Four-particle contribution to the conformal dimension in the RT-model. The second 
column shows the difference between the values of the conformal dimension of the twist field cor- 
responding to central charges C4 = 7/10 and c 2 = 1/2. The third column shows the numerically 
computed four-particle contribution to the conformal dimension for #0 = 20. 



the two particle contribution can be easily calculated to 
~, . 2cos(^) f°° 



2cos (£) r MM (l+^oELiCosh^e f0 £ cosh( ^sinh(^)sinh^ 



2 

"12 



(Eticosh(^)) 2 coshf^ 



COS \n) 

(6.83) 

where fo = m-ro is a dimensionless parameter proportional to the mass scale. From a physical 
point of view, we expect this contribution to produce a function with a plateau at exactly 
A^O) = { n ~ ~)i which is the value corresponding to two copies of the Ising model or 
c = 1. The four particle contribution is also quite simple to compute, as only few form factors 
contribute. This is because, for each copy of the model, the only non-vanishing four particle 
form factors of the energy-momentum tensor are F^l^ (^l) ^2; #3 ; #4) and all other form factors 
that can be obtained from this one by changing the particle ordering. This form factor was 
given explicitly in [48] . Together with our solution (|5.7ip and the ansatz (|5.52p it gives the four 
particle contribution 

4 V ° ; 256w 3 e<V- ./_«, Q^ =1 cosh(^)) 2 

^sh (V) cosh ) (cosh (fta) - cos (=)) (cosh (ft*) - cos (=)) 

(6.84) 

where /^i,/U 2 = + and ^3,^4 = — ■ This contribution, when added to (|6.83|) should bring the 
value of A 7 " closer to the expected one, which is obtained by setting c = 6/5 in (|3.19p . Our 
numerical results shown in Fig. 1 clearly demonstrate this to be the case for various values of n. 
As t — > —00 the functions A(t) all approach the expected value (|3.19p for n = 2,4, 6 or 8 and 
c = 6/5 with great accuracy. 
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n=2 



n=4 






Figure 1: The function A(t) := A 7 "(t) with t = 21og(f) and f = mr. In these figures we show 
the behaviour of A 7 "(t) along the renormalization group flow, from the infrared to the ultraviolet 
fixed point, for different values of the resonance parameter a. Our results are consistent with 
(|3.19p and c = 1 for the first plateau and (|3.19p with c = 6/5 for the second plateau. 



In Fig. 1 we also see that the function A(i) exhibits two finite plateaux along the renormal- 
ization group flow, which in numerical terms exactly correspond to the two particle and four 
particle contributions. The position at which the second plateau emerges changes as a function 
of a, as is also illustrated in the figure. An entirely similar behaviour was found in [37] for the 
c-function of the same model and in [51] for its effective central charge. A detailed physical 
interpretation has been given there. 

Unfortunately, the errors on the 6 particles contribution were too large to give acceptable 
results. 



7 Conclusions 

This work has been inspired by the relationship between branch point twist fields and the ground 
state entanglement entropy of 1+1 dimensional integrable QFTs. This connection provides one 
of the main motivations to study the properties of this kind of twist field. Particularly valuable 
information is provided by the form factors. They can be directly employed to generate a low 
energy (infrared) form factor expansion of the Renyi entropy (jl.4p or, as we have seen here, to 
extract the ultraviolet conformal dimension of the twist field. From a mathematical viewpoint 
one can also consider the equations (|2.8|) - ()2.12p in their own right, investigate the properties of 
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their solutions and try to find patterns as the particle numbers are increased. 

In this paper we have constructed higher particle solutions to (|2.8p - (|2.12p for two particular 
models and have tested those solutions against the A-sum rule. Our computations have revealed 
a number of interesting features: first, although the solution procedure and equations have many 
similarities with those for other local fields, it is considerably harder to find higher particle 
solutions for the twist field. This is mainly due to the increased number of poles the form 
factors have within the extended physical sheet. As a consequence, even for simple models 
it does not seem possible to find the nice closed determinant formulae found for example in 

[521231 HS1 0S| . 

Second, this extra difficulty becomes particularly clear for the Ising model for which one 
would expect to be able to find more general results. In fact, the Ising model case suggests that 
solving equations (|2.8p - (|2.12p may not be the most effective way to construct the form factors 
of twist fields. In line with this, we would like to investigate whether or not other approaches, 
such as angular quantization [65j ESI EZ] may be more suitable. In particular, the expression 
(|2.15p found in [23] was obtained in a very natural way by the latter method. 

Finally, for the RT-model we have noted that solutions to the form factors equations for 
branch point twist fields are generally not unique. This lack of uniqueness is not unexpected. 
This is because this geometric picture of the twist field as an object that connects the various 
sheets in a Riemann surface is not the only feature that characterizes the twist field. As its name 
indicates it is mainly characterized by a branch cut. One may also change the features of the 
twist field by putting other fields at the corresponding branch point. The expectation is that 
such changes will give rise to other twist fields with higher ultraviolet conformal dimensions. Our 
analysis of the RT-model, including the investigation of the cluster decomposition property of 
form factors, has revealed that some of these other twist fields correspond to non primary fields 
at conformal level and are likely to be related to composite fields involving the entropy-related 
twist field and other fields of the theory. We have found that for the RT-model and generally 
any model with a single particle spectrum, the most general solution for the 2fe-particle form 
factor of the twist field depends on k free parameters. It would be very interesting to count 
the number of twist fields systematically by counting the number of solutions to (|2.8p - (j2.12p . 
Also, we would like to investigate the conformal dimensions of all these extra fields and generally 
identify their counterparts at conformal level. 

Concerning the numerical computations performed here, our aim has been to test the few 
form factor solutions obtained for two theories: the roaming trajectories model and the SU (3)2- 
homogeneous sine-Gordon model. Both share the appearance of staircase patterns for the associ- 
ated effective central charges [381EI]- For the HSG-model the same pattern has been reproduced 
for Zamolodchikov's c-function [55\ [37] and for the conformal dimensions of certain local fields 
[37] . Our numerics demonstrate that such pattern is again reproduced for the conformal dimen- 
sion of the twist field which exhibits two plateaux at A 7 " = ^ i n ~ n) an< ^ = 5o ( n ~ n)~ 
For the RT-model we focused on the first and second steps in the staircase pattern only, corre- 
sponding to A T = (n — ^) and A 7 " = ^fg (n — ^), respectively. 

An interesting conclusion that can be drawn from our numerics, specially for the SU (3)2- 
HSG model, is that the function A^(ro) given by (|3.17p behaves exactly as 

where c(tq) is Zamolodchikov's c-function. Should this identity be exact, it would mean that 
the function A^(ro) is positive definite and monotonically decreasing (as a function of t = 
21og(mro)), just as c(ro). It is however not obvious why (|3,17p should have these features. 
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Clearly, they tell us something fundamental about the nature of the correlation function (0(r)7~(O)) 
and the branch point twist field. It would be very interesting to investigate this further, par- 
ticularly its implications (if any) for the bi-partite entanglement entropy of integrable QFTs 



Acknowledgments: The authors are grateful to Patrick Dorey, Benjamin Doyon and Andreas 
Fring for their feedback on this manuscript. 



A Explicit formulae for Q^i, X2, £3, £4) and K±(x\, X2, £3, x±) 

The constants in ()4.39p are given by 

/ .7t. ,(B-4)tt. .(2 + B)tt. 

7 = 2(l + 2cos(-))sec(-)sin(^-^)sin(^-^), (A.86) 

S = -(4 C os(-)-cos(-)-cos(^-^)jsec(-), (A.87) 

,7T% , 7T , , (1? — 1) 7T, 

= 3 + 2cos(- +4 cos — )cos( v n y ), A.88 
re 2n 2n 

f ,7i\ ,2vr N ,f 3 ,(£-1)71-, ,(5 - lWA 

£ = 2 l + 3cos (-)-cos — +8cos — cos(- — + cos(- — ) , (A.89 

\ re re 2n 2re n J 

( ,7I\ ,27T. ,37T. / ,7T . \ ,(5— lW. 

A = -2 6 + 6cos(-)+4cos( — ) + cos( — ) + 1 + 2 cos(-) cos( v 



n n rev re / re 

-4 (5eos(-) + 2cos(-) +cos(-)j cos(L__^), (A .90) 

,7I\2 / /7 T N .27T, , 7T 3 ,(5-1)71, 

p = 8cos(-) 3 + 3cos(-) + cos( — ) + 8cos( — ) cos(- — ) 

n \ n n 2n 2re 

,7I\ 2 / /7T.\ ,{B — lW. 

+8cos - l + 2cos(- cos- — ). (A.91 

re V re / n 

All the constants above are real for B real and they remain real when B = 1 — , as one would 
expect. 

The constants in (|4.40p are given by 

A = (A.92) 

(l + 2cos(^)) 3 



2(l + cos(|)) 
(l + 2cos(£)) J 



5 = 7^— ^ ( A -93) 



2(2 + cos(^)) 

C = i luZ^L a.94 

(l + 2cos(^)) 2 ' 

16cos 4 m 

D = - ML -s l (A.95) 



(1 + 2oob(=)) 



^ = 8cos(^) 2 (3 + 6cos(^) + cos(f)) 
(l + 2cos(^)) 3 
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F 



G 



2(2co S (^) + co S (gf)) 
(1 + 2 cos©) 3 
16cos(^) 2 cos(^) (2 + cos©) 



(A.98) 



(A.97) 



H 



I 



(l + 2cos©) 2 
128cos(^) 4 cos© 2 

(1 + 2 cos©) 3 
8 (cos(^) + cos(t)) 2 (3 + 2 cos© + cos(f )) 



(A.100) 



(A.99) 



(l + 2cos©) 3 



J 



256cos(© 4 cos© 4 
(l + 2cos©) 3 



(A.101) 
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